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Abstract 

Optical antennas link objects to light. Here, we analyze metal nanorod antennas as cavities with 
variable reflection coefficients to derive the interaction of dipolar transitions with radiation through 
the antenna modes. The presented analytical model accurately describes the complete emission 
process, and is summarized in a phase-matching equation. We show how antenna modes evolve as 
they become increasingly more bound, i.e. plasmonic. The results illustrate why efficient antennas 
should not be too plasmonic, and how subradiant even modes can evolve into weakly-interacting 
dark modes. Our description is valid for the interaction of nanorods with light in general, and is 
thus widely applicable. 
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Optical antennas improve the interaction of an object with optical radiation by means 
of a near-field coupling. The object absorbs and emits light through the antenna modes 

m 

. Metallic nano-particles are especially suited as optical antennas because they support 
confined plasmon modes that respond strongly to light [s], Q]. With optical antennas, the 
electronic transitions of quantum emitters, such as molecules and quantum dots, can be 
controlled. Excitation and emission rates are enhanced [s], g], Q], the spectral dependence 
shaped 8], and the angular emission directed 

0,1 

To understand optical antennas, and how they differ from conventional antennas, the Mie 

n 

solutions are available for ellipsoids [8|, and extensive numerical studies are performed for 



other shapes 



Si 



Perot) cavities 



More intuitively, antennas have been described as resonators or (Fabry- 
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If the wave vector in/along the cavity is known, 
[iqI . However, the functionality 



the position of the resonant modes can be determined 
of an antenna is not given by the value of the resonance length or wavelength alone, but by 
how its modes interact with a local object and with radiation. 

In this letter, we derive the interaction of dipolar transitions with radiation through 
optical antenna modes by treating the antenna as a cavity resonator. The obtained analytical 
model accurately describes the emission characteristics: the radiative decay rate, quantum 
efficiency and angular emission. We use the model to study the continuous evolution of the 
antenna modes from perfectly-conducting antenna theory to quasi-static plasmonics. 

We consider an elongated antenna of physical length Lp with a central section, of constant 
cross-sectional shape and size, that supports a charge density wave with wave vector: 



k = k' + ik". (1) 

The antenna is terminated at both ends, forming a resonator, which we model as a two- 
mirror cavity. Fig. [^a). The model developed applies to any cross-sectional shape, assuming 
k is known. 

The waves originate from a source at position a along the antenna axis. Figure [1] shows 
three different sources: an electric dipole, a magnetic dipole and a transmission line. The 
dipoles represent electronic dipolar transitions; transition rates are proportional to the emit- 
ted power. Electric and magnetic dipoles differ in the direction of the induced waves on 
opposite sides of the source. The transition rate for an electric (magnetic) dipole depends 
on the electric (magnetic) mode density, i.e. the magnitude of the impedance. The trans- 

2 




i z=a k 
I — • > 



(d) 







FIG. 1: (a) The antenna (total length Lp) is a rod of constant cross-sectional shape and size that 
supports a charge-density wave (wave vector k), and is terminated at both ends. It is driven by a 
local source at position z = a. We model the antenna as a ID cavity with length L and amplitude 
reflection coefficient r. In a comparison with numerical calculations, we study cylindrical antennas 
with radius R and hemispherical ends, (b)-(d) Three local sources and the direction of the waves 
induced: (b) electric and (c) magnetic dipole, (d) transmission line 

mission line resembles a magnetic dipole, but the fraction of the energy fed into the antenna 
is determined by impedance matching instead {20 . 
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22l |. As a result, dipolar transitions 



dominantly excite different modes than standard center-fed antennas. 

Resonant modes are expected for ph ysic al antenna lengths that are shifted from the 
multiples of vr/Zc' by a constant value jr.llGl. 13]. When modeling the antenna as a cavity, 



this displacement can be introduced by a positive phase shift upon reflection 1 14 . [iGj or 
by an extended cavity length [lol, [2^ . The two corrections give the same resonant length, 
but are otherwise not equivalent. We choose to set an extended length L = Lp + Lc and a 
real- valued reflection coefficient r. 

To derive the resultant current distribution I{z, a) we do not distinguish between conduc- 
tion and displacement currents and assume a one-dimensional (ID) sinusoidal distribution. 
A superposition in complex notation for time-harmonic waves gives, for ^ z < a, 



and, for a < z ^ L, 



liz, a) 



liz, a) 
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The initial amplitude of the induced wave Iq depends on the type of dipole, its oscillator 
strength, and the three-dimensional (3D) configuration. 

The + signs in equations [2] and [3] are for electric dipoles, the — signs for magnetic dipoles; 
electric and magnetic dipoles couple effectively to the antenna modes at different positions, a 
result of the geometrical argument in figured! The magnetic mode density maxima coincide 
with the electric mode density minima. 

The diffracted far-field observed at ro is given by: 

Eg = ir]o^ sine / I{z,a)e-'^°"'°'^dz, (4) 



4vrro Jq 



in which tjq is the impedance and /cq fhe wave vector for the surrounding medium. The other 
components of the electric field are zero. After evaluating the integral, equation H] becomes: 



Eg 



UqEo 



re 
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in which, k\\ = k^ cos{6) is the projection of k^ along the antenna axis, and Eq = 
zT^o^oe*'""^" sine/(47rro). A = e*^" ± re"*'"'e^*''^ and B = re*'"' ± e"*''" contain the depen- 
dence on the dipole position. 

The angular emission in equation O is one of the main results of this letter. It gives a 
complete description of the interaction of the antenna with a dipole and with radiation. It 
describes the emission of the dipole through the antenna mode and, by reciprocity [9|, its 
excitation by radiation. 

Next, we use the derived model to study the main characteristics of optical antennas in 
a set of concrete examples, and compare the results to numerical simulations. We show 
how these characteristics evolve as the antenna modes become increasingly more bound, i.e. 
plasmonic. As a measure of how bound antenna modes are, we define an effective index K: 

K = k'/ko. (6) 



As a concrete case we choose cylindrical gold antennas with hemispherical ends, 
[U^a). We study three radii, R = 20, 10 and 5 nm, which give different values for K 18 
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We vary the antenna length L for a constant wavelength. As a source, we choose an electric 
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Lk'/Tl (=j) 

FIG. 2: Evolution of the radiation resistance Rj.ad{L) for increasingly bound antennas, i.e. increas- 
ing K. The optical antennas {K = 1.7,2.9 and 5.9) are intermediate cases between the limits of 
perfect electrical conductors {K = 1) and quasistatics {K = 100 >> 1). Resonant modes occur if 
Lk' /tt = j, with j an integer. Lines: ID model. Circles: 3D Numerical calculations for cylindrical 
gold antennas in vacuum (CST MicroWave Studio), Fig. [TJ Parameters Aq = 826.6 nm. 3D 
Numerical: Eau = —29 + 2.0i. Electric dipole at 5, 2.5 and 1.25 nm from the antenna end, placed 
on and oriented along the antenna axis. ID model: Electric dipole at a = and r = 1. For K = 1 
and K = 100, k" = 0. For R = 20, 10 and 5 nm: k/ko = 1.7 + 0.048i, 2.9 + O.lli and 5.7 + 0.26i, 
and Lc = 54, 26 and 12 nm. 

dipole at the antenna end, because it effectively excites all relevant resonant modes: a = 
{Lc is added right of the dipole). 

To study the radiation damping, we define a radiation resistance for r = 1 as: 
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Rrad = '2P/Irnax^ (7) 

with P the total emitted power and Imax the maximum of (Eqs. [2] and [3]). The 

radiation resistance gives the radiation damping per unit amplitude in the resonator; it is 
independent of the total amphtude and is a characteristic of the spatial distribution of the 
mode. 

The evolution of the radiation resistance with increasingly bound modes is illustrated in 
figure [21, which shows Rrad as a function of L for the three optical antennas, together with 
the limiting cases of i^' = 1 (thin perfectly-conducting antenna), and large K (quasi-static 
limit). We make the following four observations. First, the analytical results match the 3D 
numerical calculations. Second, the modes excited by electric dipoles at a = differ from 
transmission-line center-fed antenna modes [23 [.Magnetic dipoles at a = L/2 do reproduce 



the results for center-fed perfectly-conducting [23] {K = 1) and carbon nano-tube [2^ {K = 
100) antennas. Third, unlike for K = 1, the radiation resistance for optical antennas does 
not increase with increasing length; the waves are bound. Fourth, the radiation resistance 
decreases with increasingly bound modes, i.e. increasing K. 



In the limit of >> 1, equation [5] yields Rrad oc l/K^, and since K oc 1/R [l8|, this 
also implies Rrad oc R^. The resonances are scale invariant only if Lc oc R, or equivalently if 
the reflection phase is constant, which explains previous unexpected calculation results |17] . 
and definitions jl^. 

We label the resonant modes j = 1,2... with L = jTi/k'. Even and odd modes evolve 
differently; the radiation resistance of even modes diminishes with increasing K. These 
modes have ant i- symmetric current distributions and no net dipole moment. For small 
antennas, i.e. K >> 1, these modes become subradient; the antenna modes scatter in all 
directions, opposite-oriented current elements cancel and the radiation resistance tends to 
zero. The nanorod antenna is not a simple Fabry-Perot cavity. The dependence of the 
radiation resistance on the antenna length L implies that the reflection coefficient r is not a 
constant. We relate r to the radiation resistance: 

^(^^) = ^ - Rrad/2 ^ 
Z + Rradf^ 

in which Z is the real part of the antenna wave impedance, which is taken as a free parameter 
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Lk'/7i (=j) Lk'/n (=j) 

FIG. 3: The radiative transition rate T^ad relative to the rate for j = 1 (a) and the quantum 
efficiency r] (b) for the three optical antennas. All parameters as in figure 2, but r from equation 
El with Z = 165, 250 and A50n. 

here. Equation [8] is obtained by equating the reflection loss in the cavity model with the 
antenna radiation. It assumes: that all radiative loss is due to reflection, which is partly 
justified for bound waves; that the dissipation is small so that Imax is an approximate 
measure for the current at all positions; and that Rrad depends weakly on r. 

With r defined, the radiative transition rate Trad (oc P) can be compared quantitatively 
for the different resonant modes. Fig. [3t^a). The relative values for Trad agree well with 
the numerical results. If r is taken constant instead, larger deviations are observed (e.g. 
for R = 10 nm, the error for the j = 2 peak is 8% for r{L) and 18% for r constant). 
While for K = 1 all modes are pronounced, even modes disappear with increasing K. If 
the loss is dominated by dissipation, subradiant modes evolve into dark modes with small 
Trad, despite the decreased damping due to low Rrad- By reciprocity, a small Trad implies 
low field enhancements j^; these dark modes interact weakly with radiation. For larger K, 
higher order modes are weaker compared to the j = I mode; for low Rrad the dissipative 



7 



losses {k", per length) dominate the radiative losses (r, per reflection/roundtrip), and Trad 
decays quickly with increasing L. 

The balance between radiative and dissipative (r^^) rates gives rise to a quantum effi- 
ciency rj = Trad/i^rad + Tnr)- The intrinsic efficiency of the dipole emitter is taken as unity. 
The constant additional dissipation due to the proximity of the dipole to the metal is not 
included in the model and is subtracted from the numerical results. Thus, rj is the antenna 
efficiency and sets an upper limit to the quantum efficiency of emission through the antenna 
modes. 

The efficiency generally decreases with K, particularly for the even modes, because the 
radiation resistance decreases. Fig. El^b). Clearly, an efficient antenna should not be too 
plasmonic and should in general operate away from the quasi-static small-particle plasmon 
resonance. In applications where efficient conversion into a photon is not required, large- 
K subradiant modes with low radiation damping, and thus narrow line-widths, can be 



advantageous. Examples are sensors 2J] and spasers 25 1. 



The angular emission (Eq. [5]) describes under which angles the antenna emits and can be 



effectively excited. Unlike previous OD models [26||, our model gives the emission patterns 
of higher order modes in good agreement with numerical calculations (Fig. H]). Even modes 
do not interact with radiation perpendicular to the antenna expected by symmetry 



arguments 
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271] . The lossy nature of the modes introduces asymmetry, which reveals the 
position of the dipole. Higher order modes can give multi-lobed patterns with an odd or even 
amount of maxima for odd or even modes respectively. U K >> 1 then k\\ ± k ^ ±k, and 
the 6 dependence of the denominator terms can be neglected. The emission is then a sum 
of three dipole terms: Eq, Eqc'^^^^^ and Eqc'^^^^"", with the latter contribution negligible for 
strong modes. The emission is given by two dipoles at the antenna ends, making a nanorod 



similar to a 2-slit configuration, a nd g iving a basis for the intuitive picture of scattering of 
the mode at the antenna ends [uJ, Qj 
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We summarize the interaction of the modes with radiation in a phase-matching equation 
for nano-particles: 

kii + (2m + l)kL = k'. (9) 

In which ki = tt/L, and m = 0, 1,2... approximately give the maxima of interaction (Fig. 
H]). Modes that do not give solutions for equation [9] do not interact effectively with radiation 
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FIG. 4: Angular {6) emitted power for modes j. R = 20 nm, all other parameters as in figure [3l 
ID model: line. 3D numerical: circles (blue). Phase matching, Eq. [9) dots (red). 



under any angle, and are subradiant/dark modes. Odd modes always give at least one 
solution, m = (j — l)/2. Even modes only yield solutions if Lko > vr, a condition depending 
only on L and as expected for a diffraction problem. 

To conclude, the derived model accurately describes the interaction of dipolar emitters 
with radiation through nano-rod modes. The antenna properties are primarily governed by 
a single parameter K = k' /kg that describes how plasmonic the antenna modes are, and are 
summarized in a phase-matching equation. Although here we focused on the evolution of 
the emission properties for increasingly bound waves, the model applies to all interactions 
with any spatio-temporal beam and is equally valid for field enhancement and scattering 
problems. The results are thus widely applicable and might lead to further insights and 



design rules for optica 



in/from nano-rods 
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antennas, nano-rod spasers [25|, and generally for coupling light 
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